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Abstract 

A recently developed local-density current functional formalism for con- 
fined Bose-condensed superfluids requires visco-elastic spectra which are 
defined through a finite-frequency extension of the dissipative coefficients 
entering the linearized hydrodynamic equations of the two-fluid model. We 
evaluate these spectra for a superfluid with contact interactions in the col- 
lisionless regime at finite temperature, by working to first order beyond the 
Bogolubov approximation. We find that within this approximation all the 
visco-elastic spectra take the same value aside from simple multiplicative 
factors. 
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1 Introduction 



Several experiments probing the dynamics of confined Bose- condensed alkali 
vapours have been reported in the recent literature [1-6]. In this connection, 
and following the trace set by theoretical work on the dynamics of inhomoge- 
neous electron fluids in the normal state @, f§, a Current-Density Functional 
Theory (CDFT) has been developed for an inhomogeneous superffuid in the lin- 
ear response regime |§. In particular, in the case of weak space dependences 
the CDFT equations take the form of a finite-frequency extension of the lin- 
earized hydrodynamics of the two-fluid model, taken at the local equilibrium 
densities of superfluid and normal fluid. The local-density exchange-correlation 
kernels entering the CDFT equations are complex functions of frequency describ- 
ing visco-elastic behaviour, i.e. each kernel includes both a dissipative spec- 
trum and a frequency- dependent modulus which are related to each other by a 
Kramers-Kronig relation. These two components yield in the low-frequency limit 
the dissipative friction coefficients of the linearized hydrodynamics of the super- 
fluid and the thermodynamic values of the velocities of first and second sound, 
respectively |T(J [Tl| . 

In further work [12| the bulk and shear viscosity spectra of a Bose-condensed 
fluid have been evaluated by a method previously developed for degenerate elec- 



tron fluids [13, [L4|. In this approach the equation of motion for the current 



density is treated by an approximate decoupling to account for the processes of 
excitation of correlated pairs of particles. One then recovers from the generalized 
bulk viscosity at resonance the early results of Beliaev |0| and Popov [16] for 
the width of the phonon against decay into two phonons. 

In this work we present a calculation of all the visco-elastic spectra for a dilute 
Bose-condensed fluid at finite temperature, following an approach which was de- 
veloped in early calculations of the van Hove dynamic structure factor S(k, u) by 
Wong and Gould |T7|, rj| . The interactions between the particles are described by 
means of a Fermi model potential, i.e. through a parameter v = 4na/m where a 
is the s-wave scattering length and m is the particle mass (h = 1). We express the 
spectra in terms of regular (i.e. proper and irreducible [|19[]) functions and then 



evaluate them by treating the interactions to first order beyond the Bogolubov 
approximation. As was emphasized by Wong and Gould [|I7[] , by dealing with the 
regular quantities this approach has two main advantages: (i) the single-particle 
spectrum and the density fluctuation spectrum coincide, since all the elementary 
excitations are derived from the zeroes of a suitably defined dielectric function 
(see also [pD[ ); and (ii) particle conservation is ensured by treating directly the 
longitudinal- current response, to which the density-density and density-current 
response functions are related via the continuity equation. Therefore, one is 
working in an approximation which is both gapless and conserving. 

In Section 2 we point out that within this so-called one-loop approximation 
the five visco-elastic spectra of the superfluid are determined by the values of 
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four proper "building blocks" . These are a condensate kernel related to the shift 
in local chemical potential due to interactions in the fluid away from equilibrium, 
two noncondensate kernels describing the regular part of the longitudinal and 
transverse current- current response, and a cross condensate-noncondensate vertex 
function. In Section 3 we evaluate these functions and find that the values of the 
visco-elastic spectra differ only by simple multiplicative factors. We also recover 
the equation-of-motion results ]l2j for the viscosity spectra. A brief discussion 
concludes the paper in Section 4. 



2 Building blocks for the visco-elastic spectra 

We use the standard definition Xab( t — r ',t) = — iO(t) ([A(r, t), B(r', 0)]) for a 
causal response function relating to two operators in the Heisenberg representa- 
tion, with 9(f) the Heaviside step function and (...) the expectation value in the 
equilibrium ensamble. We also introduce the Fourier transform XAB(k,u). The 
visco-elastic spectra of the superfluid are as follows 0: 

9^3 M = lim-p-Sxv.v, (*,£*>) (1) 
KCiM = S&M = }™ ^Qft&.ftw) (2) 



3? 

and 



4 

C2H + 7^{u) 



—m 2 uj^ r 



2 

^M=Jim-^*x£(M) • (4) 

Here, v s is the superfluid velocity and j is the current density, which has been de- 
composed into its longitudinal (L) and transverse (T) components. The notation 
chosen for the spectra serves to remind us that in the limit uo — > they reduce 
to the dissipative coefficients entering the linearized hydrodynamic equations of 
the superfluid (see e.g. Khalatnikov 



2.1 Vertex function, regular current response and self- 
energies 



Wong and Gould JT7| studied in detail the longitudinal current-current response 



function in the RHS of Eq. (y), which is related to the van Hove structure factor 
by 

2 

% x f.(k,Lu) = -^ r S(k,u;) . (5) 

However, they also treated the single-particle Green's function, which has the 
meaning of an amplitude- amplitude response for the Bose fluid, and introduced 
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vertex functions relating to it the amplitude-density and amplitude-current re- 
sponse functions. Evidently, the vertex functions connect the condensate with 
the noncondensate (see also 0). As we shall see below, their treatment allows 
the evaluation of all the spectra in Eqs. ©-(H). 
The single-particle Green's function is defined by 

G^ u (k,uj) = (K M ;aL))^ , ( 6 ) 

where = for fi = +, a kM = a^ k for \i — — and we use the notation 

/■oo 

((A;B)) u = -i dtexp(iut)([A(t),B(Q)]) . (7) 
Jo 

We similarly define the current-amplitude response function 

C^^^dJ^al^ (8) 
and the longitudinal and transverse current-current reponse functions 

X%(h,u) = ((k-J k ;k-J k i)) u (9) 

and 

The response functions involving the superfluid velocity are 

X VsVs {k,uj) = (k 2 /Am 2 n )f3^G^(k,u)(5 v (11) 

and 

Xy,.(k,u) = (k/2m^r )C;(k,u)^ (12) 

where /3 M = sgn(/x) and no is the condensate density. Summation over repeated 
indices is implied in these expressions. 

By virtue of the fact that only the long-wavelength limit of the response func- 
tions is needed in Eqs. ([[])- ((4|), their evaluation reduces to that of their proper 
equivalents. An explicit proof is given in Appendix A within the one-loop ap- 
proximation. We denote proper quantities by a tilde and regular (i.e. proper and 
irreducible) quantities by a (r) superscript. Following Wong and Gould jl7f , we 
define the vertex function A^ by the relation 

= KjGvh ■ (13) 
The diagrammatic analysis that they give leads to the results 

*jj = xT + K^^K (14) 

and 

Xl = xJ M • (15) 



Furthermore, the proper Green's function is expressed in terms of the proper 
self-energies M^ v through the Dyson equation as 



where 
and 



Gp, = Np,/D , (16) 
N a f3 = 6 a p[sgn(-a)u + e k - + sgn(a)sgn(/5)M_ a _ /3 (17) 



D = (uj-A) 2 -(e k - f i + S-M 2 )(e k - f i + S + M 2 ) . (18) 

In the denominator we have used the symmetrized combinations of the self- 
energies: S = (M++ + M )/2, A = (M ++ - M__)/2 and M 2 = M+- = M_ + . 

Finally, e k = k 2 /2m and fi is the chemical potential. 

The evaluation of the visco-elastic spectra thus involves the quantities that 
we have introduced in Eqs. (^l)-(ID- We proceed to treat these quantities in a 
perturbative framework. 



2.2 One-loop approximation 



In the Bogolubov approximation the excitation spectrum reduces to a single un- 
damped mode at each wavenumber, with a dispersion relation changing continu- 
ously from phonons at long wavelengths to free particles at high momenta. The 
condensate density n is equal to the total density n and the total current co- 
incides with the superfluid (or condensate) current. The values taken by the 
quantities that we have introduced in Section 2.1 are = fi^ = vn, M 2 (0) = 0, 
i(o) = o, A*(°) = ^k(3 fl /2m and x L ' T(r)(0) = 0. 

Broad visco-elastic spectra arise in the so-called one-loop approximation, in 
which one treats the corrections to the Bogolubov results at first order in the 
dimensionless coupling-strength parameter g = ((47ra) 3 n) 1//2 . From Eqs. (|13|)- 
fll8|) we calculate to this order of approximation the proper functions that are 
needed for the evaluation of the visco-elastic spectra, with the following results: 



(i) 



/ >2 a 2 



(5<i) + m« _ ^) 



(i) 



(19) 



V 



(1) 



(20) 



~L(1) 
v •• 



k 



+ (S® + M 2 (1) - ^) + — ^ (2uAM + k 2 (S^ - ^ 



(21) 



5 



and 



Y- (1) 
A JJ 



(22) 

In these equations 5^ is the vector (1, 1) and n' = n — hq is the depletion of the 
condensate. We have used the convention of summation over repeated indices 
and expressed wavenumbers in units of ^/Arrna and energies in units of 4irna/m. 
In these units the coupling strength parameter is g = Ana = 1/n. 

It is evident that the evaluation of Eqs. (|i~9l) -(f22|) involves seven quantities, 
that is (S^ + M^-fi^), A£(%, Af >/? M , x L / )(l) and xJS r){1) . 

The relevant one-loop diagrams for these functions are shown in Figure 1, where 
the single-particle propagators represent Bogolubov quasiparticles. In the present 
units each condensate line carries a factor g^ 1 ^ 2 and each current vertex carries 
a factor g 1//2 . The self-energy diagrams in Figure 1 include two-quasiparticle 
excitation processes via scattering against the condensate, in addition to the usual 
Hartree and Fock terms as for a normal fluid. The vertex function accounts for 
damping of longitudinal-current fluctuations via excitation of quasiparticles and 
scattering against the condensate. Finally, the attenuation of both longitudinal 
and transverse current fluctuations proceeds via two-quasiparticle processes. 

We show in Appendix B that the quantities (S^ 1 ' — fx^) and K^ 1 '^ do 
not contribute to the visco-elastic spectra in the long wavelength limit. There- 
fore, only the four quantities (S« + M 2 (1) - A*<%, x-/ ){1) and xj (r)(1) 
remain to be calculated in that limit. The first of these is a condensate ker- 
nel contributing to all three longitudinal spectra, while the second is a cross 
condensate-noncondensate kernel involved in the longitudinal current spectrum 
and in the correlations between longitudinal current and superfluid velocity fluc- 
tuations. Finally, two noncondensate kernels enter the attenuation of longitudinal 
and transverse current fluctuations. 



3 Calculation of visco-elastic spectra 

The calculation of the diagrams shown in Figure 1 yields the following expressions 



for the building blocks of the visco-elastic spectra at zero temperature |T7], [18 

d 3 p 



+ M 1] - = -7T / ^A p A p+k ^ - „° +k - u° p 
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r d 3 p 
(2tt) 3 


3(xf )(1) ) 




r d 3 p 
(2vr) 3 


and 






ex./ ~T(r)(lK 


= -- 1 
16 J 


, d 3 p 
(2tt) 3 



1 



A n A 



p A p+k 



0J r 



p 



(p.kf 



(Ap — A p+k ) 
ApA p+ k 



(23) 



p-k + -k)(X p - A p+k ) 6(u> - <, k - u° p ) , (24) 



(25) 



-8{u - - W J) . (26) 



6 



In these equations we have defined u>® = yp 2 +p 4 /4 and A p = e p /u; p , being 
the Bogolubov quasiparticle frequency in reduced units. 

The thermal factor (l+/ p+ k+/p), with f p = (exp(/?u; p ) — l) -1 and/? = 1/ksT, 
must be inserted in the integrands in Eqs. (|23|)- (p6|) at finite temperature, on 
account of detailed balance in the process of decay of an excitation of energy u 
into two quasiparticles. In addition, there appear other contributions associated 
with the scattering of the excitation against a thermally excited quasiparticle. 
However, the corresponding integrals contain a factor (/ p +k — fp), which vanishes 
in the long wavelength limit. 

Bearing in mind the above observations, we easily find from Eqs. (P3"D-(f2l)|) 
the following expressions for the four building blocks in the long wavelength limit 
at finite temperature: 



fc-o v ^ 2 ^ ' cosh(/3c<j/2) - 1 16ttcj'(1 +p 2 /4) 



smh{(3u/2) kp 3 



and 



lim SWA^/J 1 — 

k-a 1 " ^ cosh(/^/2) - 1 48vr^(l + p 2 /4) 3 / 2 

sinh(/?cu/2) k 2 p 2 



(27) 
(28) 



LU°=iu/2 



lim Sfv-- 

fc^o VA JJ cosh(/5w/2) - 1 1607rj;(1 + p 2 /4f 

~T(r)(l) 



w°=w/2 



3hmQ(jft r > w ). (29) 



Here, uj' = du^/dp and p is defined by the condition = u>/2. 

Using these results in Eqs. (|T§|)-([2"2"D and hence in Eqs. ([TT|)-(]T3D, we obtain 



K fo(w)) = i» + ^)) = ^nmR (£(«)) = ^n 2 m 2 K (Cs(w)) 

sinh(/3u;/2) n 2 v 2 p 6 



cosh(/?a;/2) — 1 307rm 2 u; 3 co> 



(30) 



We remark that the behaviour of these spectra at low frequency is proportional 
to u 3 at T = and to uj 2 T at finite temperature. As already noticed in the 
Introduction, the result for the shear and bulk viscosity spectra agrees with that 
obtained in [0 by an equation-of-motion approach in an approximation which 
takes account of the excitation processes of two quasiparticles. 

Figure 2 shows the behaviour of the shear viscosity spectrum at zero tem- 
perature and at finite values of the temperature. We remark from Eq. (|30D that 
the spectra involve the Fourier transform v of the interatomic potential, so that 
their high-frequency behaviour depends on this potential at high momenta and 
hence on the details of the short-range interatomic collisions. Difficulties can be 



7 



anticipated in assessing this asymptotic behaviour for the Bose-condensed alkali 
vapours: in this case a realistic description of the short-range potential (a) is not 
Fourier transformable, and (b) gives a bound state which is not to be included 
in the study of the metastable condensate state. 

A simple way to regularize the contact potential that we have assumed in the 
present calculations would be through the introduction of a high-momentum cut- 
off jl^]. For a reasonable value of the cut-off the results shown over the frequency 
range illustrated in Figure 2, as well as the simple proportionality of the various 
spectra which is seen in Eq. (]30f) , would remain approximately valid. 

The Hilbert transform of the visco-elastic spectra gives the frequency- dependent 
corrections to the thermodynamic values of the moduli determining the first and 
second sound velocity. Evidently, regularization of the interaction potential will 
be necessary before evaluating these corrections. 

4 Concluding remarks 

We have calculated the visco-elastic spectra for a dilute Bose-condensed fluid 
in the collisionless regime at finite temperature. The motivation for this work 
has been provided by the theory developed in ||, showing that these spectra as 
determined for a homogeneous Bose fluid allow a local density treatment of the 
dynamics of a weakly inhomogeneous superfluid. 

We have assumed contact interactions between the particles and treated them 
to first order beyond the Bogolubov approximation. Within this scheme we have 
found that in the long-wavelength limit the spectra are determined by four proper 
building blocks, to which we have attributed meaning in terms of dynamic corre- 
lations between condensate and noncondensate. Explicit evaluation of these four 
functions has shown that within this approach the visco-elastic spectra are in fact 
determined by a single function of frequency. However, the high-frequency part 
of the spectra and the visco-elastic shifts in the sound velocity moduli will be sen- 
sitive to details of short-range interatomic collisions beyond the contact-potential 
scheme. 

We have briefly indicated how our calculation of the shear and bulk viscosity 
spectra for the superfluid is related to an equation-of-motion treatment of the 
dynamics of current fluctuations. The latter approach has also been applied to an 
electron fluid in the normal state. Since we have treated the dissipation spectra in 
the superfluid as arising from the interactions between Bogolubov quasiparticles, 
we have found that the current-density fluctuations can only decay into pairs of 
longitudinal excitations. In contrast, decay channels invoving transverse current 
fluctuations also contribute to damping in a normal fluid. 

Experimental studies of higher frequency excitations in the collisionless regime 
would be useful for quantitative tests of the theory. The presently available exper- 
iments on confined Bose-Einstein condensates refer to low-frequency excitations, 
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in a dynamical regime where collisions begin to become important depending on 
the thermodynamic state of the fluid. The theory of the transition to a collisional 
regime, where the dynamical behaviour of the fluid is governed by a Boltzmann 
equation at frequencies below an average collision rate, remains to be explored. 
A start in this direction was made in the early work of Popov ||16|| . 
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Appendix A. Reduction to proper quantities 

In this Appendix we demonstrate that within the one-loop approximation the 
response functions introduced in Eqs. (|9|)-(|T2|) reduce to the corresponding proper 
response functions in the long wavelength limit. 

We start from the superfluid velocity response in Eq. (|11]). From the dia- 
grammatic analysis given by Wong and Gould J[7| for the single-particle Green's 
function one easily finds 



e 0) _, v ( k ~ 2 

e De \lj M oj 



OJ 



A) 



(31) 



where e has the meaning of a dielectric function in terms of the proper density- 
density response function and is its regular part. A lengthy but straightfor- 
ward calculation yields the one-loop result 



OJ 



-k c> 



OJ 



OJ 



f)2 



+ 



2e k 
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4* a*(D 5 

^-ojff* 

(i) 



+ 4(o; 2 -4)xf ) -fc 2 xf )(1) + L/ 



uu* - ei 



+ 



4el 



oj^-uf) 



-9 



/ 2 2\~L{1) . 



2uj 



jj 



OJ 



(32) 



in terms of the functions entering Eqs. (|I^)-(^T|). The long wavelength behaviour 
of these functions, from Eqs. (p3|)-(p6|) and from the results given in Appendix B, 
is as follows: all the self-energy terms are of order k°, the vertex-function terms 
are of order k 1 and the current response functions are of order k 2 . Hence, from 
Eq. (H) it follows that 



lim ^(P^G^Pt, 

k— >0 



lim ^(P^G^P^ 



di) 



(33) 
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Similarly for the J — v s response function we have 



k. 



k. 



u 2 (u 2 - e\ 



-3 



L.2~L(r)(l) 



- 2fcAf)^ + 2(5 (1) + M 2 (1) -/i«; 
and one immediately finds 

lim9f(^)« = lnn9f(C^)W. 



fe^O 



(34) 



(35) 



Finally, for the current response function one can use directly Eqs. (3.31) and 

(36) 



(3.32) in the paper of Wong and Gould [|17] to obtain 



1- CV L{1) 

lim ^Y:: 



lim 



A similar property is evidently valid for xjj- 



Appendix B. Some further one-loop results 



We complete the calculation of the one-loop spectral functions in supplement to 
Eqs. (E3b-(EB). We have 



3i« 



S(S-/i) (1) 
and 



3(AW (1) 



d 3 p 
(2^ 

d 3 p 



(27T) 



ApAp +k + - 1 



P- k + -k 



A 



p+k 



1 6(u — uj 



i P\ 
p+k - ^ P J 



(37) 



5{u - u° p+k - u° p ) (38) 



(39) 



From these expressions it is easily seen that in the long wavelength limit the 
quantities ( |3"7| ) and ( Pffi ) are of order k° and the quantity ( |3"9"D is of order k 1 . 
Therefore, their contributions to Eqs. (0)-(^l|) vanish in this limit. 
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Figure captions 



Figure 1. One- loop diagrams for the proper self-energies M^, vertex function 
and current-current response functions Xjj- The full lines are Bogolubov 
single-particle propagators, while the wiggly lines represent the condensate and 
the triangles denote current vertices. The interaction lines are shown in dashes. 

Figure 2. Viscosity spectrum 77(0;) at zero (solid line) and finite temperature 
(dashed T = 3a;o/8, dot-dashed T = 3cjo/4) as a function of frequency in units 
of Uq = 4:7van/m. 
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